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Abstract 



Fractional calculus [jl], |^ is a branch of classical mathematics, which deals with the 
generalization of operations of differentiation and integration to fractional order. Such 
a generalization is not merely a mathematical curiosity but has found applications 
in various fields of physical sciences. Various differential equations of physics, such 
as wave equations, diffusion equations [|, ||, H, etc, have been generalized to include 
fractional orders. Recently there has been a surge of activity involving the use of 
fractional calculus in studying phenomena relating to fractal structures and processes. 

Fractals P, 0] are sets having fractional dimension. They appear naturally in 
nonlinear and nonequilibrium phenomena in different forms and context. Frequently 
attractors of chaotic dynamical systems are fractals (e.g. fractal dust or nowhere dif- 
ferentiable fractal curve). We note that the graph of a nowhere differentiable functions 
is generally a fractal set. Also typical paths of a Brownian particle or Feynman paths 
in Quantum Mechanics p|, |^ are nowhere differentiable curves. There are processes, 
e.g. Levy process, which generate a fractal set. Processes which have evolution taking 
place at time instances lying on fractal sets can be used to model various transport 
phenomena, e.g., transport in spongy or disordered materials, transport in chaotic 
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Hamiltonian systems. 

Ordinary calculus is inadequate to characterize and handle such irregular curves 
and surfaces. Also we felt a need to develop calculus which will describe phenomena 
(such as diffusion) on fractals. It is an interesting question to ask: whether the 
formalism of fractional calculus can provide such a framework. The aim of the thesis 
is to explore this possibility. And our answer is in the affirmative. 

The first chapter of the thesis contains motivation for the work that we have 
carried out. It also gives the overview and defines the scope of the thesis. 

The second chapter gives the introduction and brief review of the related fields. 
Which includes a survey of the relevant part of the fractional calculus. Definitions and 
examples of fractals are also given. A classical notion of fractional differentiability 
has also been briefly reviewed. 

The third chapter starts with the deflnition of the local fractional derivative(LFD), 
for order g, < g < 1, along with its motivation. The LFD of order g of a function 
f{y) is deflned by 

iDV(y) = lim^^iP^^M o<g<l (0.1) 

where the derivative on the RHS is the Riemann-Liouville fractional derivative [jl], 0, 
viz., for < g < 1 

d'fix) 1 d f- f{y) 



[d{x - a)]" r(l - g) dxJa (x - y^^^' ^'^''^^ 

This deflnition is then generalized to include all positive orders g. Further we have also 
generalized this deflnition for functions of several variables. All the deflnitions have 



been illustrated with the help of examples. After this introduction of the definitions 
a fractional Taylor expansion has been derived which involves these LFDs. Various 
consequences of this local fractional Taylor expansion have been discussed. 

The fourth chapter deals with the question of fractional differentiability properties 
of continuous but nowhere differentiable functions. It starts with the introduction to 
nowhere differentiable functions. Then a specific historically important example of 
Weierstrass function is considered. Using the definition of LFD introduced in the last 
chapter we have proved the striking result that this function is locally fractionally 
differentiable upto the critical order 2 — s where s is the box dimension of the graph 
of the function. This is followed by some general results which establish a connection 
between the critical order and the local Holder exponent of the function. 

Use of LFD in studying regularity properties of the irregular functions has been 
considered in the fifth chapter. The chapter begins with the introduction and exam- 
ples of multifractal functions. This follows by the demonstration, with the help of 
examples, of the fact that LFD can be used to detect isolated singularities as well as 
singularities which are masked by stronger ones. Further it has been shown that the 
LFD can be used even in the case where singularities are not isolated and local Holder 
exponent is not same at all points. We have used a specific example of a selfsimilar 
multifractal function constructed by Jaffard for this purpose. 

The sixth chapter introduces a new kind of equation called local fractional differ- 
ential equation (LFDE) which is an equation involving LFDs. This is a generalization 
of usual differential equations to fractional order keeping the local nature intact. A 
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simple example of such an equation, viz., 

E>lf{x)^g{x), (0.3) 

is considered. It is found that the solution to this equation does not exist when 
g{x) — Const.. But it admits a finite and nontrivial solution if g{x) has fractal 
support. It is argued that these LFDEs are suitable to describe phenomenon taking 
place in fractal space and time. As one possible application of LFDE to physical 
situation a local fractional Fokker-Planck equation has been derived starting from 
first principles. A specific case of this equation is solved and found to give rise to a 
subdiffusive solution. 

The seventh and final chapter concludes the thesis. Along with the concluding 
remarks we discuss the further question which have come up. 
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Chapter 1 



Motivation and Overview 



"One might encounter instances where using a function without derivative 
would he simpler than using one that can he differentiated. When this hap- 
pens, the mathematical study of irregular continua will prove its practical 
value. " 

. . . Perrin 

Till the second half of the last century, it was believed that a continuous function 
must be different iable at least at some point. This view changed in 1872, when 
Weierstrass constructed a function which was continuous but nowhere differentiable. 
Subsequently, a host of instances of other such functions was developed (historical 
comments are summarized in For about three decades, these functions were 

considered as pathological cases and were used only as counter examples. As can be 
seen from the above quotation (see [^), it was Perrin who recognized their relevance 
to the physical world. In his study on Brownian motion, he realized that the path of 
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a Brownian particle can be modeled by such irregular functions. 

Subsequent advances in nonlinear and nonequlibrium phenomena revealed many 
more examples of the occurrence of such functions in physical theories. These in- 
clude, among others, various generalizations of the Brownian motion [0, ^, Q (e.g., 
fractional Brownian motion), Feynman paths in quantum mechanics, attractors of 
some dynamical systems ||^, |T0|, velocity field of turbulent fluid ||^, |^, etc. Before 
proceeding to the discussions of these irregular functions we would like to point out 
the developments in the mathematical theory of sets of fractional dimension showed 
that generally the graphs of such functions are fractals (hence we may alternatively 
call these functions as fractal functions). 

• Jean Perrin (who won the Nobel Prize for his work on Brownian motion) pointed 
out that the path followed by a Brownian particle is irregular on all length scales 
and hence, one can not draw a tangent to this path. It was shown later that 
the path followed by the Brownian particle in with embedding dimension 
ii^ > 2 is a fractal with dimension 2. Another similar example comes from 
quantum mechanics where typical Feynman paths are differentiable nowhere. 
Observed path of a particle in quantum mechanics was shown to be a frac- 
tal curve with Hausdorff dimension 2. The work on Brownian motion (see 0]) 
showed that the graphs of components of Brownian paths are nowhere differen- 
tiable and have dimension 3/2. Brownian motion was later generalized |^, Q 
to fractional Brownian motion. This process is known to give rise to graphs 
which have dimensions between 1 and 2. 
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Turbulent fluid system is another source of irregular functions. In these systems, 
the inherent nonlinear nature in the turbulent regime is believed to give rise 
to the irregularity. Passive scalars advected by a turbulent fluid have been 
shown 1^ to have isoscalar surfaces which are highly irregular, in the limit 
of diffusion constant going to zero. Velocity fields of a turbulent fluid at low 
viscosity p[ have been shown to be multifractal. 

In the study of dynamical systems theory, attractors of some systems are found 
to be continuous but nowhere differentiable M. We will consider one specific 



example [0, ITO[ of the dynamical system which gives rise to such an attractor. 



Consider the following map. 

Xn+l = 2Xn + yn mod 1, 

Vn+l = Xn + Vn mod 1, (1.1) 
Zn+1 = XZn + COs(27rX„) 

where x and y are taken mod 1 and z can be any real number. In order to 
keep z bounded, A is chosen between and 1. The equations for x and y are 
independent of z. The x — y dynamics are chaotic and are unaffected by the 
value of z. It was shown in that if A > 2/(3 + y^), then the attractor of this 
map is a nowhere differentiable torus (i.e., it has same topological form as a 
torus but is nowhere differentiable) . A simpler example of a map whose repeller 
is a graph of a nowhere differentiable function is considered in O. 



A variety of natural and industrial processes lead to formation of rough surfaces. 



Crystal growth, vapor deposition, grinding, erosion, fracture are a few examples 
of processes leading to such surfaces. To understand the dynamics of growth 
and formation of rough interfaces is one of the challenging problems. Formalism 



based on concepts of fractals [0, [T^, |T3| has become a standard tool in the 
study of such surfaces. 

Many problems of physical interest involve solving a partial differential equation 
with Dirichlet boundary conditions. A variety of methods exist to solve such 
equations with smooth boundaries. But solving such equations with nowhere 
differentiable surfaces as a boundary presents a great difficulty. Such problems 
arise when one wants, for instance, to find field near a charged rough surface |jT^ 



or to study vibrations of a drum with fractal boundary [15 



The world is full of such nonlinear and nonequilibrium phenomena. However, 
the complicated nature of these phenomena make them analytically intractable. One 
of the reasons for the lack of availability of analytical methods can be traced to 
almost ubiquitous appearance of fractals in these phenomena. It is clear from the 
examples given above that the graph of the fractal functions can arise as attractors 
or repellers of dynamical systems, as paths of particles performing random motion, 
etc. Moreover, general fractal sets appear at numerous places in nonlinear theories. 
There are many interesting problems in which one has to deal with the phenomena 
taking place in fractal space or time. For instance, one may consider the diffusion on 
fractals, where underlying space is a fractal. Such a model is used to describe diffusion 
in disordered systems. Diffusion in which transitions take place at time instances 
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which he on a fractal set, can be used to model diffusion in presence of traps. The 
usual differential calculus, which is used to describe phenomena in Euclidean space 
and time, is inadequate to handle these cases. 

It is interesting to investigate whether fractional calculus, which generalizes the 
operation of derivation and integration to fractional order, can provide a possible 
calculus to deal with fractals. In fact there has been a surge of activity in recent 
times which supports this point of view. This possible connection between fractals and 
fractional calculus gives rise to various interesting questions. For example, one may 
ask that since an n-fold integral gives an n-dimensional volume, whether a suitably 
defined 'a- fold' integral {a noninteger) would give rise to a-dimensional measure. Or, 
how would one generalize a formula for finding length of a smooth curve to one which 
finds fractal measure of a continuous but nowhere differentiable curve? These lead to 
a more specific problem of studying fractional differentiability properties of nowhere 
differentiable functions, and, the possible relation of the order of differentiability with 
the dimension of the graph of the function. 

The aim of the thesis was to investigate the issues discussed above. The strategy 
we employed was as follows. We first concentrated on the fractal (irregular and 
possibly nowhere differentiable) functions and study their fractional differentiability 
properties. Our attention was focused on the local scaling property of these functions. 
This led us to construct a quantity called local fractional derivative. It was then used 
to study pointwise behavior of irregular functions. Finally, we demonstrated that 
these local fractional derivatives can be used to construct equations to deal with 
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phenomena involving fractals. It can also be seen as a first step towards development 
of much sought after calculus for describing phenomena taking place in fractal space 
and time. 

Multifractal measure has been the object of many investigations H, |T6|, O, |18 



19] , |20| , |2T| , [2^ , with a host of interesting applications. Its importance also stems 



from the fact that such measures are natural measures to be used in the analyses 
of many physical phenomena |12, 1^. Rigorous mathematical results concerning 



multifractals were developed in [^. However, it may so happen that the object one 
wants to understand is a function (e.g. a fractal or multifractal signal), rather than 
a set or a measure. For instance, one would like to characterize the velocity field of 
fully developed turbulence. Studies on fluid turbulence have shown the existence of 
multifractality in velocity fields of a turbulent fluid at low viscosity [||. 

In the light of these needs there has been a renewed interest, chiefly among math- 
ematicians, in studying pointwise behavior of multifractal functions. New tools have 
been developed to study local behavior of functions, which will be of practical use. 
Notable among them is the use of Wavelet transforms ["M, 25, 26, 27 1 which has some 
success in this effort. In this thesis we develop and use local fractional derivatives 
to study pointwise behavior of multifractal functions and point out some advantages 
over the existing methods. 

The further plan of the thesis is to begin by surveying some relevant fields in the 
next chapter. With this introduction, we proceed in the thesis with the characteri- 
zation of the local scaling property of the fractal functions. In order to do this, we 
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modify the definition of the usual fractional derivative and introduce the notion of lo- 
cal fractional derivative (chapter ^). It is shown, in the same chapter, that these local 
fractional derivatives appear naturally in a fractional Taylor expansion. In chapter ^ 
we study the local fractional differentiability properties of the nowhere differentiable 
functions. Choosing the Weierstrass function as a prototype example, we show that 
this function is locally fractionally differentiable at every point, upto some critical or- 
der which lies between and 1. This critical order is related to the fractal dimension of 
the graph of the function. We then demonstrate (chapter ^) that the local fractional 
derivatives can be used to study pointwise behavior of the multifractal functions. In 
chapter ^ we introduce the local fractional differential equations. These equations 
are expected to describe phenomena in fractal space-time. The algorithm for the 
solutions of special cases of these equations has also been discussed in this chapter. 
As an application, we derive the fractional analog of the Fokker-Planck equation. A 
special case of this equation is then solved and shown to have a subdiffusive solution. 
The thesis ends with the concluding remarks. 
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Chapter 2 
Introduction 



The purpose of this chapter is to present a brief introduction to some of the back- 
ground material needed for the theme of the thesis. Since it deals with applications 
of fractional calculus to fractals, introduction to both these fields is essential. In the 
first section we consider the fractional calculus, giving various definitions of fractional 
derivatives and integrals, their properties and few examples. In the next section we 
deal with fractals. We list desirable properties of the definition of dimension and, 
in particular, discuss two definitions in detail. We also discuss some applications of 
fractals in the physical systems. In a separate section we review earlier applications of 
fractional calculus to scaling phenomena. The last section of this chapter is devoted 
to a brief review of classical notions of fractional differentiability. 
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2.1 Fractional Calculus 



Thus it follows that d^x will he equal to x\J dx : x, an apparent paradox, 



from which one day useful consequences will be drawn. 



. . . Leibniz 



Fractional calculus |l|, |^ is a branch of classical mathematics which deals with 



generalization of the operations of derivation and integration to fractional order (even 
to irrational and complex orders - in this sense it is a misnomer). Though little 
known, the subject is by no means new. Historically the interest in the subject dates 
all the way back to the year 1695, when Leibniz made the above prophetic remark 
in a letter to L'Hospital, quoted in [|[] (this was the era during which Newton and 
Leibniz independently developed the ideas of classical calculus). During the next two 
centuries the subject slowly matured through the works of the mathematicians Euler, 
Lagrange, Laplace, Abel, Liouville, Riemann, to name but a few. 

2.1.1 Definitions 

There exists a multitude of definitions of fractional derivatives and integrals, which 
have different origins and are not necessarily equivalent. We shall list few of them here. 
Let's start with the definition given by Grunwald (see 0]), which directly extends and 
unifies the notions of quotient of differences and Riemann sums. According to this 
definition a differintegral of order g of a function f : M ^ M is given by 



[d{x - a)]i 



d'f 



lim 

N-*oo 




N-1 



j=0 



E 



r(j - q) 
r(j + 1) 



f{x-j 



x — a 




(2.1) 
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where q, a & IR and T{p) is the usual gamma function defined by the integral 

T{p) = / xP-^e-'^dx, p>0, (2.2) 





and by the relation r(p + 1) = pr(p). For g > the above formula yields derivatives 
of order q, and for g < it gives integrals of order q. We shall follow this convention 
throughout. It may be noted that the fractional derivative depends on the lower limit 
'a '. The above definition involves the fewest restrictions on the functions to which 
it applies and, unlike some other definitions introduced below, does not use ordinary 
derivatives and integrals of the function. However, in practice it is somewhat difficult 
to use this definition, except in the cases of simple functions. 

The most frequently used definition of a fractional integral is the Riemann- 
Liouville definition. According to this definition, a fractional integral of order g of a 
function / is given by 

d'f(x) 1 r fiv) 



[d{x-a)]i T{-q)Ja (x - 
where the lower limit a is some real number. The fractional derivative of order g is 

d'^fix) 1 d^ r fiy) 



r 7 ,-, dy for n-l<q<n (2.4) 



[d{x — a)]i V{n — q)dx'^Ja {x — yY 
_ {X - a)-^^VW(a) 

" h r(-g + + 1) + [d{x - a)Y-- ^ l<q<n. (2.5) 

As is clear, this definition uses the concepts of ordinary derivatives and integration. 
Since it amounts to evaluating an integral, it is more convenient to use. It can 
be shown that the Grunwald definition and the Riemann-Liouville definition are 
equivalent to each other. 



Hermann Weyl defined fractional integral of order g of a function / by 

This definition is suitable for periodic functions as it leaves the periodicity of the 
function unaffected, unlike in the cases of previous two definitions. Because of this 
property, Zygmund has used this definition extensively in working with trigono- 
metric series. Many more definitions of fractional derivatives and integrals can be 
found in references 0]. 

2.1.2 Properties 

In this subsection we review some of the properties of the operations of differentiation 
and integration to arbitrary order |^, |^ . 

Linearity and homogeneity 

From the definitions introduced above it is clear that the differintegral satisfies 

[d{x - a)]i [d{x - a)]? [d{x - a)Y ' ' ' 

and 



[d{x-a)]i [d{x-a)]'i' 
Change of scale 

When the argument of the function is scaled by a factor f3, the differintegral satisfies 

[dx]'} ^ [di/3x)]i' ^ ' 



19 



For a more general formula, with nonzero lower limit a, see [Q. 



Leibniz's rule 

Osier proved a general result for a differintegral of a product of two functions. It 
is given by 

dVg] _ ^ T{q + 1) dmfd-r+jg 

where 7 is arbitrary. He used 0, H this result to generate certain infinite series 
expansions which interrelate special functions of mathematical physics. 

Chain rule 

If is an analytic function and / a sufficiently differentiable function, then the chain 
rule for fractional differentiation is given by ||T| 

where J2 extends over all combinations of nonnegative integer values of Pi, P2, . . . , Pn 
such that J2k=i = n and J2k=i = This is a generalization of Faa de Bruno's 
chain rule for integer order derivatives. The complicated nature of this formula makes 
it less useful. 

Differintegration term by term 

Rules for distributing fractional derivative and integral operators over an infinite series 
are the same as those for ordinary derivatives and integrals. Let us first consider the 
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case of g < 0. Suppose the infinite series of functions J2 fj converges uniformly in 
< |x - a| < X, then 11 



M oo OO Mf 

and the right-hand side also converges uniformly in < |a; — a| < X. For g > 0, if 
the infinite series J2 fj as well as the series J2 d.'^fj/[d{x — a)]'' converges uniformly in 
< Ix — a| < X, then 

oo oo jq f 

T.f. = T.T17^\T,, ^>0, (2.13) 



[d{x - a)]" f^^ [d{x - a)]" 

for < |x — a| < X. 

Fractional Taylor series 

In [0 Osier generalized the Taylor series to include fractional derivatives. He proved 
a very general result. We consider a special case which is called the Taylor- Riemann 
series: 

where h ^ zq and f{z) is analytic. Osier used this generalization to study some 
special functions of mathematical physics. Please note that this series contains terms 
with negative powers of {z — zq) and hence it is not suitable for use in standard 
approximating schemes. 
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Composition rule 



Here we study the relationship between two operators 



r^^^-.,11.+Q - (2.15) 



[d{x - a)]'? [d{x - a)]^ [d{x - a)] 

When / = 0, 

d'i ci«[0] rf^+«[0] 



(2.16) 



[d{x - a)]'? [d{x - a)]Q [d{x - a)]^+Q 
and the composition rule is trivially satisfied. When / 7^ 0, the composition rule is 
satisfied if 

d-Q d^f _ 

■'~ [d{x-a)]-Q [d{x-a)]Q ~° (^■^'^^ 

and if this condition is not satisfied then we have 

d'^ d^f d'i+^f 



[d{x - a)\i [d{x - a)]^ [d{x - a)] 



q+Q 



[d{x - a)Y+Q - a)]"Q [rf(a; - a)]Q 

It can be shown that when Q < 0, the composition rule is satisfied for an arbitrary 
differintegrable function /. In fact if f{x)/{x — a)'^ is finite at a; = a, then the 
composition rule is valid for Q < p + 1. 

2.1.3 Examples 

In general, finding fractional derivative or integral of a function can turn out to be 
very complicated. In this section we consider few examples of familiar functions for 
which the job is relatively easy 0, |^. 
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The function 

First we find the fractional integral of the function using the Riemann-Liouville 
definition, i.e., we have, with a = 0, 



d^xP 
dxi 



yP 



T{-q)Jo (a;-y)9+i 



dy for g < 0. 



(2.19) 



With a change of variable, we get 



d^xP 
dxi 



X 



r(-g) 



■p-q 



r(-g) 



/ uPil-u^-^du, for g<0 
Jo 

B{p + l,-q), q<0,p>-l, 



(2.20) 
(2.21) 



where B{p, q) is the beta function. Therefore fractional integral of a x^ is, for p > — 1, 

d^xP T{p + 1) ^_ 
— x^ 



dx'i V{p — q + 1) 



q < 0, p> -1. 



(2.22) 



To find fractional derivative we use equation ( p.4|) and obtain 



d'^xP 
dxi 



dx'^ 



X' 



p—q+n „i 



V{n — q) Jo 



uPil - uY~^-^du 



for < g < n 



r(p + i) 

T{p-q + l) 



X- 



p-i 



q>0,p> -1. 



(2.23) 



The following limiting behavior of the abovestated formula may be noted, viz.. 



lim 



d'^xP 



+0+ dxi 



ii q<p or q = p + n 

T{p +1) if g = p 
oo otherwise. 



(2.24) 



This observation will be of crucial importance later in the thesis. 
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The exponential function exp(x) 
We use the power series expansion 



expfx) = y — -, (2.25) 



,t^or(j + i)' 

which is valid for all x, and differintegrate it term by term. Therefore we get 



ci^expfa;) _„ ^ x^ 



dx" fr'o T{j -q + 1) 

7*(-g,x) 



xi 

where 7* (a, b) is an incomplete gamma function 

The logarithmic function ln(x) 

We use the Riemann-Liouville integral and get 
d'^\n{x) 1 /"^ \n{y) 



(2.26) 



dxi T{-q)Jo (x-y)''+i 

a;~^ln(x) dv x^'^ f'^ \n{l — v)dv 



,9+1 



(2.27) 



7(— g) Jo li—Q) Jo 

where we have made use of the substitution v = {x — y)/x. Evaluating the above 
integrals we get 

d'^ ln(x) x"^ 



dx'^ 7(1 — g) 
where ipi^) is defined by 



[ln(a;) - 7 - ^/'(l - q)] (2.28) 



Using equation (|2.4| ) this result can be extended to g > 0. 
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Trigonometric functions 

The differintegration of sine and cosine functions results in complex hypergeometric 
functions and is given in The purpose of this section is to note some relations, 
among fractional derivatives of sine and cosine functions. 

If we choose f{x) = sinx and a = in the equation ( |2.5[ ) we get 

d'' sin{x) d'^~^ cos{x) 

dx'i ~ dxi-^ ' ^ ' ' 

where < g < 1. Similarly, when we choose f{x) = cosx and a = we get, for 
< g < 1, 

d'^ cos{x) x~'^ d''~^ sm{x) 



dx'i r(-g + 1) dxi~^ 
These relations will be of use later in our discussion. 



(2.31) 



2.2 Fractals 



2.2.1 Introduction 



Loosely speaking, fractals are sets or objects having non-integer dimensions. Enor 
mous body of literature already exists on this very active field |2^, ||, pT 



Conventionally, the dimension of a object is usually a non-negative integer (though 
there are situations where negative dimensions are also sensible |2^), and it is defined 
through the number of coordinates needed for complete specification of the given ob- 
ject. Therefore, when one wonders about the possibility of existence of entities which 
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have non-integer dimension, one has to look for a different definition of dimension 
which does not depend on the coordinates. 

The intuition for such a definition comes from the way one measures length, area or 
volume in practice. To make the basic idea of the definitions of fractional dimensions 
clear we present this procedure explicitly. In order to find the length of a smooth 
curve, one takes a measuring stick of length say S and makes tick marks at regular 
intervals on the curve with the help of this stick. Then the length L of the curve is 
the product of N{S), the number of such marks, and S. In order to make an finer 
estimate of the length, one makes 5 smaller. One expects that as S approaches a limit 
the quantity L remains bounded and approaches a limit. In similar manner in order 
to find area (volume) of a set it is covered with squares (cubes) of sides 5 and the 
number N{S) of such squares (cubes) needed to cover a given set is counted. Then the 
area (volume) is defined as A — lim5_^o 5'^N(5) {V — hm5_^o 5^N{5)). Notice that the 
power of the 6 above corresponds to the dimension of the sets. Notice also that for a 
planar piece, the length is infinite, the area is finite and the volume is zero. The value 
of the dimension where the content measure of the set jumps from infinity to zero is 
the characteristic dimension of the set. This intuition will provide the foundation in 
the following formal definitions of fractional dimension. 

2.2.2 Definitions 

There exist numerous definitions of fractional dimension, each having its own advan- 
tages and disadvantages. We shall consider two important definitions of dimension 
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below. We shall follow the treatment of ref. p4|. Before going to the definitions let 



us list some desirable properties, as discussed in which a definition of dimension 
should have. Let F C iR". 

1. Monotonicity: If E G F then dimE < dimF. 

This is obvious since if we add some points to a given set its dimension should 
not decrease. 

2. Stability: dim{E [j F) = max{dimE , dimF) 

3. Countable stability: dim(lJi^i Fi) = snp^^^^^ dimFj 

The above two properties imply that if we combine a set with other sets having 
lower dimensions, the dimension of the resultant set does not change. 

4. Geometric invariance: If iR" is transformed with a transformation such as a 
translation, rotation, similarity or affinity, the dimension of the set F C iR" 
should not change. 

5. Lipschitz invariance: If f : F M"^ is a bi-lipschitz transformation, i.e. 

ci\x-y\ < \ f{x) - fiy)\ < C2\x-y\, x,y e F, 
where < ci < C2 < oo, then dim/(F) = dimF. 

6. Countable sets: If F is finite or countable then dimF = 0. 

This is again expected since the dimension of the singleton set is zero and from 
the countable stability any countable set should have dimension zero. 
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7. Open sets: If F is an open subset of IR" then dimF = n. 

8. Smooth manifolds: If F is a smooth m-dimensional manifold then dimF = m. 
Of course, not every definition of the dimension satisfies all the above properties. 

Hausdorff dimension 

Mathematically this is one of the most important definition of the dimension, since 
it is defined for any set. A major disadvantage of this definition is that it is not 
algorithmic. 

Let us first define some quantities which are needed in the definition of the Haus- 
dorff dimension. If [/ is a non-empty subset of iR", the diameter of U is defined as 
\U\ = sup{|a; — y\ : x,y E U}. If {Ui} is a countable collection of sets of diameter at 
most 6 that cover the given set F, i.e. F C [JiZi Ui, we say that {Ui} is a 6-cover of 
F. 

Now we define a quantity, for any 5 > 0, 



where s is a non-negative number. In the above equation we consider all covers of 
F by sets of diameter at most 6 and take the infimum of the sum of the sth powers 



reduced. Therefore, l-Ll{F) is a monotonic function, and approaches a limit (finite or 
-|-oo) as 5 — * 0. Now we define the Hausdorff measure by 




(2.32) 



of the diameters. As 5 decreases, the number of permissible covers of F in ( p.32| ) is 



W{F) = \\m^Hl{F). 



(2.33) 
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This limit exists for any subset F of iR". 

It can be proved that if T-C'^{F) is finite for s = Si then T-C'^{F) = for s > Si. 



This imphes that there exists a critical value of s at which 7i^(F) jumps from oo to 
0. This critical value is the Hausdorff dimension of F. 



It is a customary p3| to normalize the Hausdorff measure by multiplying it with 



_ |r(i/2)l' 

' r(i + s/2)' ^ ' 

which is a s-dimensional volume of a ball of unit radius. 
Box dimension 

This is the most popular definition of dimension since it is easy to use both mathe- 
matically and numerically. 

In order to define the box dimension consider the collection of cubes in the 5- 
coordinate mesh and Ns{F) be the smallest number of such cubes required to cover a 
non-empty bounded set. Then the lower and upper box dimensions of F are respec- 
tively defined as 

dimBF = lim,^o^^n^ (2.35) 
-logd 

\ogNs{F) 



dime-F = lim^^o ; (2.36) 

-logd 



The box dimension of F is defined as 



dimeF = dimpF = dimeF (2.37) 
whenever the two limits are equal. 
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This definition of dimension satisfies most of tlie aforementioned desirable prop- 
erties of dimension, except that it is not countably stable and hence countable sets 
may not have box dimension 0. For example, consider a set of rational numbers in the 
interval [0, 1]. If one covers this set by a partition with interval 5 then the N{5) = 1/5 
and hence the box dimension of the set is 1. But the set of rational numbers is a 
countable set, i.e., it is a countable union of sets of dimension zero. Hence one would 
expect the set of rationals to have dimension 0. In spite of this paradoxical result 
this definition is widely used by the physicists because of its simplicity, under the 
assumption that such cases are pathological in physical world. 

It may be noted that in the definition of the box dimension we had essentially 
covered the set with boxes of fixed size whereas in the Hausdorff dimension we allow all 
sizes smaller than some 5. This is the crucial difference between the two definitions. 
We now see that because of this difference the Hausdorff dimension of the set of 
rationals in [0, 1] is 0. Let us cover the set of rationals in [0, 1] in following manner. 
Since this is a countable set we can label each rational number by a positive integer 
say k. Now cover A;th rational number by an interval of length 5/2^ . Then the sum 
(5^/2'^* is bounded by K5^ (for some positive constant K) which goes to 0, for s > 0, 
in the limit 5^0. Hence the quantity T-C^{F) is for any s > 0. 

2.2.3 Examples 

We consider here two simple and well-known examples of fractal sets, i.e., the Cantor 
set and the von Koch curve. Graph of the Weierstrass function discussed in chapter ^ 
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provides still another example. Many more examples can similarly be constructed. 
They are found, for example, in [^, |^. 



The middle-third Cantor set 

The middle-third Cantor set is constructed from a unit interval [0, 1] by a sequence 
of deletion operations. Let Eq be the interval [0, 1]. Let Ei be the set obtained by 



deleting the middle-third of Eq (see figure |2Tl| ) . Therefore the set Ei consists of two 



intervals [0, 1/3] and [2/3, 1]. The set E2 consisting of four intervals [0, 1/9], [2/9, 1/3], 
[2/7,7/9], [8/9, 1] is obtained by deleting middle thirds of the above two intervals. 
We can continue in this way, deleting the middle thirds of each interval in Ek-i to 
obtain Ek, ad infinitum. The limiting set we get is called the middle third cantor set. 
Notice that the set E^. consists of 2^ intervals of length 3"^^. That means we need 2^ 
intervals of length 3"^^ to cover the set at kth level. Therefore the box dimension of 
the set is clearly log 2/ log 3. 

The von Koch curve 

In order to construct the von Koch curve we again start with the unit interval [0, 1]. 
Now we remove the middle third segment of this interval and instead of throwing 
it off we replace it, as shown in figure p.2| , by the other two sides of the equilateral 
triangle based on the removed segment. In this way we get the set Ei which consists 
of four segments. By applying the same procedure to each segment of Ei we obtain 
the set E2, and so on. The limiting set we get is called the von Koch curve. It may be 
noticed that the set E^ consists of 4^^ segments of length 3^^^. Therefore the dimension 
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E2 



F 



Figure 2.1: The middle third Cantor set 

of the von Koch curve is log 4/ log 3. 

These two examples of fractals, though simple, display many typical characteristics 
of fractals. The set F in both these examples is self-similar, i.e., it contains copies of 
itself on all length scales. The set F contains intricate detailed structure at arbitrarily 
small scales. In spite of this, the set is easy to construct using a recursive procedure. 
More importantly, the geometry of F is not easily described in standard notion, i.e., 
it is not the locus of points that satisfy some simple geometric condition, or it is not 
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a solution of any simple equation. 



2.2.4 Occurrences of Fractals 

The notion of fractals has found innumerable applications in pure as well as in applied 
sciences. The present thesis itself contains lot of such applications. We have already 
mentioned some examples of fractals in the beginning of the thesis. In the following 
we discuss some additional physical applications. Many more applications may be 



found in references |10, 11, 12, 13 



Kinetic Aggregation 

The diffusion-limited aggregation (DLA) model |13|, Q is a simple model of a fractal 



generated by diffusion of particles in the following manner. One starts with a 'seed' 
particle which is fixed at a point. A second particle starts on the border of a circle 
around this point and performs a random walk. When it comes near the seed it 
sticks to it and a cluster of two particles is formed. Then a third particle performs a 
similar random walk and forms a cluster of three particles. This procedure is repeated 
many times, thereby giving rise to a large cluster having the shape of a fractal. DLA 
models many phenomena in nature, such as, viscous fingering, dielectric breakdown, 
snowflake growth, etc. 

Biological systems 

There is abundance of fractal structures and processes in biological systems. Human 
lungs, branching of trees, root system in plants, etc, have a selfsimilar branching 
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structure which is typical of fractals [15]. The abovementioned DLA model has also 
been used to understand the shape of a neuron, growth of bacterial colonies, etc. 
Long range correlations, which are believed to give rise to fractal geometry, are found 
in DNA sequences, human heart beats etc. 

Percolation 

Consider a large square lattice, where each site is occupied randomly with probability 
p and empty with probability 1 — p. At low values of p, the occupied sites form small 
clusters. As one increases p there exists a threshold Pc at which a large (macroscopic) 
cluster appears and which connects opposite edges of the lattice. This is called the 
percolation cluster. This percolation cluster is selfsimilar in nature and hence has 
fractal structure. 



2.3 Applications to the study of scaling in physical 
systems 



Recent work |]T^ ^ shows that the fractional calculus formalism is useful 

in dealing with scaling phenomena. The purpose of this section is to present a brief 
review of a few developments in this area. 
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Fractional Brownian Motion 

Mandelbrot and Van Ness have used fractional integrals to formulate fractal 
processes such as fractional Brownian motion. The fractional Brownian motion is 
described by the probability distribution Buit) defined by, 

Buit) - 5^(0) = ^ f Kit - t')dB(t'), 0<H<1, (2.38) 

where B{t) is an ordinary Gaussian random process with zero mean and unit variance, 
and with K{t — t') such that 

(t - t')"'^'^ 0<t' <t 

K{t -t') = { (2.39) 
{(t_t')^^-i/2_(_t)i/-i/2| <o. 



The use of the fractional integral kernel is apparent. Recently Sebastian has 
given a path integral representation for fractional Brownian motion whose measure 
contains fractional derivatives of paths. 

Fractional equations for a class of Levy type probability den- 
sities 



Levy flights and Levy walks [31, ^ have found applications in various branches 



of physics ||3^ ranging from fluid dynamics |3^ to polymers |3^. In T. F. Non 



nenmacher considered a class of normalized one-sided Levy type probability densities, 
which provides a distribution for the lengths of jumps of the random walker, given by 



fix) = -—x-^"'^ exp(-a/x) a > 0, x > 0. (2.40) 
1 [n) 
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It is clear that for large x, 

f{x) ~ x->'-^ /i > 0, X > 0, (2.41) 



where /i is the Levy index. In |T^] it was shown that these Levy-type probability 



densities satisfy a fractional integral equation of the form 



= (2.42) 



or equivalently, a fractional differential equation given by 

\ ' = a-^^{x^^f{x)). (2.43) 

Here, n=lifO</i<l and n = 2ifl<yU<2, etc. An interesting observation 
is that the Levy index is related to the order of the fractional integral or differential 
operator. 

Fractional equations for physical processes 

Modifications of standard equations governing physical processes such as diffusion 



equation, wave equation and Fokker-Planck equation have been suggested |20, Ell 36 



37, 38 1, which incorporate fractional derivatives with respect to time. In refs. |20 



21| a fractional diffusion equation has been proposed for the diffusion on fractals. 
Asymptotic solution of this equation coincides with the result obtained numerically. 
Fogedby |39| has considered a generalization of the Fokker-Planck equation involving 



addition of a fractional gradient operator (defined as the Fourier transform of — fc^) 
to the usual Fokker-Planck equation and performed a renormalization group analysis. 
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Transport in chaotic Hamiltonian systems 

Recently Zaslavsky showed that chaotic Hamiltonian dynamics, under certain 
circumstances, can be described by a fractional generalization of the Fokker-Planck- 
Kolmogorov equation, which is defined by two fractional critical exponents {a, (3) 
respectively for the space and time derivatives of the distribution function. With 
certain assumptions, the following equation was derived for the transition probabilities 

p{x,ty. 

d^P(x t) 9" 1 9^° 

where 

, , A'(x;At) 
A(x) = lim ' 2.45 

and 

B{x) = lim (2.46) 

with A' and B' being a*^ and 2a^^ moments respectively. The exponents a and jS 
have been related to anomalous transport exponent. 

Some more instances 



Hilfer |jT6|, [T^, |1^] has generalized the Eherenfest's classification of phase tran- 
sition by analytically continuing the order which is equivalent to using the 
fractional derivative. 
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Some recent papers [19, 20, 21, 22] indicate a connection between fractional 
calculus and fractal structure |23, |2l] or fractal processes p6|, ETl p8|. 



W. G. Glockle and T. F. Nonnenmacher |^ have formulated fractional dif- 
ferential equations for some relaxation processes which are essentially fractal 



time 28 processes. 



We note that the most of the applications reviewed in this section deal with asymp- 



totic scaling only. As noted in ^1|] , in the context of Schrodinger equation containing 
fractional operators, these fractional differential equations are approximations as they 
do not completely take into account the local scaling property of the fractals. Also, in 
some of the examples above the dynamics is governed by non-local equations owing 
to the fact that the fractional time derivatives used are non-local. This point has 



also been emphasized in |^ . In the later part of the thesis, we develop a new notion 
which takes care of these problems. 



2.4 Weyl fractional differentiability 



The purpose of this section is to review briefly, for the sake of completeness, the 



classic works of Stein, Zygmund and others [y, 43, 44, 45 , who use Weyl fractional 
calculus in the analyses of irregular functions. 
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2.4.1 Holder exponent 

Let f{x) be defined in a closed interval / of real line, and let 

u{6)=u;{6:f) = snp\f{x)-f{y)\ for x,y e I \x - y\ < 6 (2.47) 

The function uj{6) is called the modulus of continuity of /. If for some a > we have 
uj{6) < C6"' with C independent of 6, then / is said to satisfy Holder (In old literature 
it is also called Lipschitz) condition of order a. These functions define a class Aq, of 
functions. Notice that with this definition of the modulus of continuity only the case 
< a < 1 is interesting, since if a > 1 for all x of J, f'{x) is zero everywhere and the 
function is constant. This definition has been generalized (see to one where the 
case a > 1 is also nontrivial. The idea is to subtract an appropriate polynomial from 
the function /. 

uj{6) = Lj{6 : f) = sup\f{x) - P{x - y)\ for x,y e I \x - y\ < 6 (2.48) 

where P is the only polynomial of smallest degree which gives the smallest order of 
magnitude for u. Recently, this definition was extensively used ^ to characterize 
the velocity field of a turbulent fluid. 

2.4.2 Summary of classical results 



Following Welland |^5| we introduce 



Definition 1 f is said to have an ath derivative, where < /c < a < k + 1 for integer 
k, if D^'^f (in the Weyl sense), with P = k + 1 — a, has k + 1 Peano derivatives at 
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Xq. That is, there exists a polynomial Pxo{t) of degree < k + 1 s. t. 

(D-^/)(xo + t) - P.0 W = t ^ 0. (2.49) 

Further if 

{- r \{D-^f){xo + t)-P,M'dtV/'' = o{p''-') P^O (l<j9<oo) (2.50) 
f is said to have an ath derivative in the sense. 



It is clear that this definition of fractional differentiability is not local. Particularly 
the behavior of function at — oo also plays a crucial role. The main results can be 
stated using this notion of differentiability and involve the classes and A*"^ which 
are given by the following definitions 



Definition 2 // there exists a polynomial Qxo{t) of degree < k s. t. /(xq + t) — 
Qxo{t) = 0{\t\°') as t ^ then f is said to satisfy the condition A^, and if 



{- r \f{xo + t)-Qxo{t)\^dty/^ = 0{p''), p^O (l<j9<oo) (2.51) 
P J-p 

f is said to satisfy the condition A^. 



Definition S f is said to satisfy the condition if for some p > 

1 fP \f{xo + t)-QxM'' 



, , -dt < oo. (2.52) 

We are now in a position to state the main results in the form of the following four 



theorems. First two of these P4| , |45[| state the condition under which the fractional 
derivative of a function exists. 
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Theorem 1 Suppose that f satisfies the condition at every point of a set E of 
positive measure. Then D°'f{x) exists almost everywhere in E if and only if f satisfies 
condition almost everywhere in E. 

Theorem 2 The necessary and sufficient condition that f satisfies the condition 
almost everywhere in a set E is that f satisfies the condition and D'^f exists in 
the sense, almost everywhere in this set. 

The following results tell us how the class of a function changes when an operation 
of fractional differentiation is performed. 

Theorem 3 Let 0<a<l, P>0 and suppose that f E A^. Then D^f^f E Aa+/3 if 
a + f3 <1. 

Theorem 4 Let < < a < 1,. Then f E A„_^ if f E ha- 

Despite their merits, these results are not really suitable and adequate to obtain 
information regarding irregular behavior of functions and Holder exponents. We 
observe that the Weyl definition involves highly nonlocal information and hence is 
somewhat unsuitable for treatment of local scaling behavior. In the next section we 
introduce a more appropriate definition. 
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Chapter 3 



Local Fractional Derivatives 



3.1 Motivation 

The definition of tlie fractional derivative was discussed in the last chapter. These 
derivatives differ in some aspects from integer order derivatives. In order to see this, 
one may note, from equation (|2.4|) , that except when g is a positive integer, the q^^ 
derivative is nonlocal as it depends on the lower limit 'a'. The same feature is also 
shown by other definitions. However, we wish to study local scaling properties and 
hence we need to modify this definition accordingly. Secondly from equation ( p. 231) 
it is clear that the fractional derivative of a constant function is not zero. Therefore 
adding a constant to a function alters the value of the fractional derivative. This is an 
undesirable property of the fractional derivatives to study fractional differentiability. 
While constructing local fractional derivative operator, we have to correct for these 
two features. This forces one to choose the lower limit as well as the additive constant 
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before hand. The most natural choices are as follows. 1) We subtract, from the 
function, the value of the function at the point where we want to study the local 
scaling property. This makes the value of the function zero at that point, canceling 
the effect of any constant term. 2) Natural choice of a lower limit will again be that 
point itself, where we intend to examine the local scaling. 

In the following section we formally introduce the concept of local fractional deriva- 
tive. 



3.2 Definition 

Definition 4 //, for a function f : [0, 1] M, the limit 

iDV(y) = lim 0<g<l (3.1) 

•'^^^ ^^y d{x-yY ^ ' 

exists and is finite, then we say that the local fractional derivative (LFD) of order q 
(denoted by ]D'^f{y)), aty, exists. 

This defines the LFD for < g < 1. It was first introduced in and later general- 
ized to include all positive values of q as follows. 

Definition 5 //, for a function f : [0, 1] M, the limit 

ID'^fiy) = hm ^ ^n-or(n+i)l ^3 2) 

^~-y [d{x-y)]'i ^ ' 

exists and is finite, where N is the largest integer for which N^^ derivative of f{x) at 
y exists and is finite, then we say that the local fractional derivative (LFD) of order 
q {N < q < N + 1), at x = y, exists. 
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We subtract the Taylor series term in the above definition for the same reason as one 
subtracts f{y) in the definition ||. We do this to supress any regular behavior that 
may mask the local singularity. 

Definition 6 The critical order a, at y, of a function f is 

a{y) = sup{g|iD^ /(?/), q' < q, exists}. 

Sometimes it is essential to distinguish between limits, and hence the critical order, 
taken from above and below. In that case we define 

"^^^^y^ = [d±{x-y)]^ • ^'-'^ 

We will assume iD^ = E)\_ unless mentioned otherwise. 

As an example consider the function f{x) = a > 0. The critical order of 
this function from above at origin is a when a is noninteger and is oo when a is an 
integer. 

The local fractional derivative that we have defined above reduces to the usual 
derivatives of integer order when g is a positive integer. In order to see this point for 



g = 1 we consider equation (|3.1|) . Since the Riemann-Liouville fractional derivative 
on RHS of equation ( |3.1|) reduces to ordinary first derivative when q = 1 (This follows 
from the relation 



d dif{x) 



X] 



dx [d{x - y)]i [d{x - 
and setting q = 0. For details, see the book by Oldham and Spanier 0, page 50. See 

also ref. 0), we get 

lD^f{y) = \im^lM:^lMl. (3.4) 
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Now since f{y) is constant its derivative is zero. Therefore after taking the hmit we 
get 

iDV(y) = (3.5) 

For N < q < N + 1 the definition of local fractional derivative is 

"^^^y^ = H K^^^iF ^ 

Now, if = n (i.e. N — n — 1 in view of the above range of q), we get 
Again since y is a constant 

^V(y) = (3.8) 

Therefore, for q — n, local fractional derivative reduces to nth order derivative. 

From this it is clear that in our construction local fractional derivatives generalize 
the usual derivatives to fractional order keeping the local nature of the derivative 
operator intact, in contrast to other definitions of fractional derivatives. The local 
nature of the operation of derivation is crucial at many places, for instance, in study- 
ing differentiable structure of complicated manifolds, studying evolution of physical 
systems locally, etc. The virtue of such a local quantity will be evident in the follow- 
ing section where we show that the local fractional derivative appears naturally in 
the fractional Taylor expansion. This will imply that the LFDs are not introduced in 
an ad hoc manner merely to satisfy the two conditions mentioned in the beginning, 
but they have their own importance. 
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3.3 Local fractional Taylor expansion 



An interesting consequence of the above definitions is that the LFDs appear naturally 
in the fractional Taylor expansion. 

Derivation of local fractional Taylor expansion 



We follow the usual procedure to derive Taylor expansion with a remainder |jTO[. In 
order to derive local fractional Taylor expansion, let 

It is clear that 

iDV(l/) = i^(l/,0;g) (3.10) 

Now, for < g < 1, 



f^'^'^-f^y^ = [d{x-y)]-'^ [d{x-y)]i^' 
1 f^~y F{y,t;q) 



T{q) Jo (x-y- t)-i+^ 
-l-[F{y,t;q) j {x - y - tf-^ dtf^-^ 



dt (3.11) 



1 r^-ydFiy,t;q)ix-y-ty^^ ^3^^^) 
Jo dt a 



m Jo 

provided the last term exists. Thus 
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i.e. 



m = m + ^^(^ - yr + ^.(^, y) (3.14) 



where Rq{x,y) is a remainder given by 



R^(^^y) = fl^^i ^^(x-y-tm (3.15) 



Equation ( p.l4| ) is a fractional Taylor expansion of f{x) involving only the lowest and 



the second leading terms. Using the general definition of LFD and following similar 
steps one arrives at the fractional Taylor expansion for < g < + 1 (provided ID'^ 
exists), given by, 

where 

R^(^^y) = fl^^l {x-y-tydt (3.17) 

We note that the local fractional derivative (not just fractional derivative) as 
defined above provides the coefficient A in the approximation of f{x) by the function 
f{y) + A{x — yY/r{q + 1), for < g < 1, in the vicinity of y. We further note 
that the terms on the RHS of eqn( p.l3D are nontrivial and finite only in the case 



q = a. Osier in ref.|]TT| has constructed fractional Taylor series using usual (not 
local in the sense above) fractional derivatives. His results are, however, applicable to 
analytic functions and cannot be used for non-differentiable scaling functions directly. 
Furthermore, Osier's formulation involves terms with negative q also and hence is not 
suitable for approximating schemes. 
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Let us consider the function f{x) = x,a > 0. Then iD"/(0) = r{a + 1) and 
using equation ( |3.16| ) at y = we get /(x) = since the remainder term turns out 
to be zero. 

Geometrical interpretation of LFD 

It is well known |^ that one can not attach any geometrical interpretation to the 
conventional fractional derivatives as one does for ordinary derivatives. Whereas, 



the local fractional Taylor expansion of section |3.3| suggests a possibility of such an 
interpretation for LFDs. In order to see this note that when q is set equal to unity in 
the equation (|3.14|) one gets the equation of the tangent. It may be recalled that all the 
curves passing through a point y and having same tangent form an equivalence class 
(which is modeled by a linear behavior). Analogously all the functions (curves) with 
the same critical order a and the same iD" will form an equivalence class modeled by 
x" (If / differs from x" by a logarithmic correction then terms on RHS of eqn( p.l3| ) 
do not make sense like in the ordinary calculus). This is how one may generalize 
the geometric interpretation of derivatives in terms of 'tangents'. This observation is 
useful when one wants to approximate an irregular function by a piecewise smooth 
(scaling) function. 

3.4 Generalization to higher dimensional functions 

The definition of the Local fractional derivative can be generalized for higher 
dimensional functions in the following manner. 
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Consider a function / : iR" IR. We define 

Hy,t) = f{y + ^rt)-f{y), vGiR", teM. (3.18) 

Then the directional-LFD of / at y of order g, < g < 1, in the direction v is given 
(provided it exists) by 

where the RHS involves the usual fractional derivative of equation ( p.4|) . The direc- 
tional LFDs along the unit vector will be called i^^ partial-LFD. 

3.5 Some remarks 

1. We would like to point out that there is a multiplicity of definitions of fractional 
derivatives. The use of Riemann-Liouville definition, and other equivalent def- 
initions such as Grunwald's definition are suitable for our purpose. The other 
definitions of fractional derivatives which do not allow control over both the 
limits, such as Wyel's definition or definition using Fourier transforms, are not 
suitable since it would not be possible to retrieve local nature of differentiabil- 
ity property which is essential for study of local behavior. Also, the important 
difference between our work and the work of p, ^ is that while we are try- 
ing to study the local scaling behavior these works apply to asymptotic scaling 
properties. 

2. It is interesting to note that the same definition of LFD can be used for negative 
values of the critical order between -1 and 0. For this range of critical orders 
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= — 1 and the sum in equation ( p.2|) is empty. As a result the expression for 
LFD becomes 

nyXy) = Urn ^ (3-20) 

An equivalence between the critical order and the Holder exponent, for positive 
values of critical order, will be proved in chapter |^. Here we would like to point 
out that the negative Holder exponents do arise in real physical situation of 



turbulent velocity field (see |12, 13| and references therein). 



3. Another way of generalizing the LFD to the values of critical order beyond 1 
would have been to write it as 

DVfa)^.im'''(^':'W-rfa)) (3.21) 

But the existence of iV*^ derivative of / at a; may not be guaranteed in general. 
Such a situation may arise in the case of multifractal functions to be treated in 
chapter 
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Chapter 4 



Fractional Differentiability of 
Nowhere Differentiable Functions 

4.1 Introduction to nowhere differentiable func- 
tions 

It is well known that a continuous function may not be differentiable at some point 
(e.g., f{x) = \x\ is not differentiable at x = 0). But, as mentioned in Chapter 
1, till the last century it was generally believed that any continuous function must 
be differentiable atleast at one point. The example of the Weierstrass continuous 
but nowhere differentiable function was a counterpoint to this belief. Weierstrass 
constructed the function 

oo 

W{t) = ^ a'' COS bH, (4.1) 

k=l 
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where < a < 1 < 6 and ab > 1 (see [T[]). While he himself and others did 
independently prove that this function is nowhere differentiable for some of these 
values of a and b, it was Hardy |l|] who for the first time showed this for all such a 
and b. 

We shall consider the following form of Weierstrass function 

oo 

Wx{t) = J2 ^^'"^^^ sin \H, A > 1 1< s < 2, (4.2) 

k=l 

where we have chosen b = X and a = A^*"^-"^. Note that Wa(0) = 0. 

The graph of this function is known to be a fractal. Fractal dimension of such 
and related curves were first studied in [Q]. The box dimension of the graph of this 
function is s [^]. The Hausdorff dimension of its graph is still unknown. The best 
known bounds are given by Mauldin and Williams |^ where they have shown that 
there is a constant c such that 



< dim// graph/ < s. 



log A 

Recently it was shown that if a random phase is added to each cosine term in (|0 



then the Hausdorff dimension of the graph of the resulting function is 2 + logo/ log 6. 
These kind of functions have been studied in detail in 0, ^, Figure ^?T| shows 



graphs of the Weierstrass function for a fixed value of A and various values of the 
dimensions. It is clear that as dimension increases the graph becomes more and more 
irregular. As the value of the dimension approaches 2 the graph tends to fill the area. 
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Weierstrass function as a limiting set of a dynamical system 

Nowhere differntiable functions arise as attractors or repellers of dynamical systems. 
One such example was cited in chapter |^. Here we discuss another example in detail. 
Consider a mapping h : —>■ IB? given by 

h{t, x) = (At, A2-^(x - sin(t))), (4.3) 

where A > 1 and < s < 2. Then we get 

h^{t,x) = h{Xt,X^-'{x-sm{t))) (4.4) 
= {X\X^-\X^-'{x-sm{t))-sm{Xt))) (4.5) 
= {XH, A2(2-^)(a; - (sin(t) + A^'^ sin(At))) (4.6) 

Therefore, continuing in the same fashion, we have the nth iterate as 

n 

h''{t,x) = (A"t,A"(2-^)(a;-^A(^-2)'^sinA^t). (4.7) 

k=l 

This shows that the graph of Weierstrass function ( [4. 21) is a repeller for h in (|4.3| ). 

4.2 Fractional differentiability of Weierstrass func- 
tion 

Now we consider the Weierstrass function as a prototype example of a continuous 
but nowhere differentiable function (which is homogeneous in some sense to be expli- 
aned later) and study its local fractional differentiability properties [0]. We begin by 
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considering the origin t = 0. We first find fractional derivative of Wx{t): 



dfi dti 

oo 



d{\HY 



k=l 



provided tlie riglit liand side converges uniformly. Using tlie relation 

d'^ sinfx) d'^"^ cosfx) 

-, 0<g<l, 



dx'i dxi 
we get 



dti d{XH)i-^ ' 



(4.8) 

From second mean value theorem it follows that the fractional integral of cos(A'^t) of 
order g — 1 is bounded uniformly for all values of X^t. This implies that the series on 
the right hand side will converge uniformly for q < 2 — s, justifying our taking the 
fractional derivative operator inside the sum. 

Also as t — s> for any k the fractional integral in the summation of equation ( |4.8| ) 
goes to zero. Therefore it is easy to see that 

iDWA(O) = lim^^^^^ = for g < 2 - s. (4.9) 

This shows that q*^ local derivative of the Weierstrass function exists and is continu- 
ous, at t = 0, for g < 2 — s. 

To check the fractional differentiability at any other point, say r, we use t' = t — r 
and W{t') = W{t' + r) - W{t) so that W{0) = 0. We have 

oo oo 

Wxit') = ^A(^-2)^sinA'=(t' + r)-5]A(^-2)^sinAV 

k=l k=l 

oo 

= 5]A("-2)'^(cosAVsinA¥ + sinAV(cosAV-l)) (4.10) 



k=l 
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Taking fractional derivative of this with respect to t' and following the same procedure 
we can show that the fractional derivative of the Weierstrass function of order q < 2—s 
exists at all points. 

For q > 2 — s, right hand side of the equation ( [4 .81 ) seems to diverge and our first 
step of taking fractional derivative operator inside the sum itself may not be justified. 
Therefor to prove that the LFD of order q > 2 — s does not exist we adopt a different 
approach. We do this by showing that there exists a sequence of points approaching 
along which the limit of the fractional derivative of order 2 — s < g < 1 does not 
exist. We use the property of the Weierstrass function 0, viz, for each t' G [0, 1] and 
< 6 < 6o there exists t such that |t — 1'| < 6 and 

cd"" <\W{t) -W{t')\, (4.11) 

where c > and a = 2 — s, provided A is sufficiently large. We consider the case of 
t' = and t > 0. Define g{t) = W{t) - ct°. 

The abovementioned property, alongwith continuity, of the Weierstrass function 
guarantees the existence of a sequence of points ti > ^2 > ••• > > ••• > such that 
— > as n — > oo and g{tn) = and g{t) > on e) for some e > 0, for all n (It is 
not ruled out that t„ may be zero for finite n). Let us define 

gn{t) = if t<tn 

= g{t) otherwise. 

Therefore we have, for < a < g < 1, 

d'^gnjt) ^ 1 d g{y) ^ 
d{t - tnY r(l - q) dtJu (t - y)i ^' 
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where tn <t < t„_i. We assume that left hand side of the above equation exists. Let 

Now h{tn) = and /i(t„ + e) > 0, for a suitable e, as the integrand is positive. Because 
of the continuity there must exist an e' > and e' < e such that h{t) is increasing on 
e'). Therefore 



Q^^^^nW I n = 1,2,3.... (4.12) 

- d{t - tnY ^ 

This implies that 

c—, -\t-t < M-j-j n = 1,2,3... 4.13 

But we know from eqn( |2.23| ) that, when < a < g < 1, the left hand side in the 
above inequality approaches infinity as t 0. This implies that the right hand side 
of the above inequality does not exist as t — > 0. This argument can be generalized 
for all nonzero t' by changing the variable t" = t — t' . This concludes the proof. 

Therefore the critical order of the Weierstrass function is 2 — s at all points. Thus 
the maximum order of differentiability is directly related to the dimensions. 



4.3 Critical order and Levy index of a Levy flight 

Schlesinger et al |^ have considered a Levy flight on a one dimensional periodic lattice 
where a particle jumps from one lattice site to other with the probability density given 
by 

Ll! — 1 °° 

P(x) = J2 ^'^{^ix, W) + 5{x, -b^)] (4.14) 
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where x is magnitude of the jump, 6 is a lattice spacing and h > uj > 1. 5{x, y) is a 
Kronecker delta. The characteristic function for P{x) is given by 

_ 1 oo 

P{k)^^^Y/;j-^^os{Vk), (4.15) 

which is nothing but the Weierstrass cosine function. For this distribution the Levy 
index is logo;/ log 6, which can be identified as critical order of P{k). 

More generally, for the Levy distribution with index the characteristic function 
is given by 

P(A;) =^ expel A;]'^. (4.16) 

Critical order of this function aXk — Q also turns out to be same as ji. Thus the Levy 
index can be identified as the critical order of the characteristic function at /c = 0. 



4.4 Equivalence between critical order and the Holder 
exponent 

Theorem 5 Let f : ^ IR he a continuous function. 

a) If 

y d\f{x) - f{y)) 

lim — ;— -T = for o < q; , 

^-^y [d{x-y)]'i 

where q,a & (0, 1), for all y e (0, 1), then dimBf{x) <2 — a. 

b ) If there exists a sequence x^^ y as n ^ oo such that 

hm — r— -^1 = ±oo for q > a , 

^ [d{xn - y)]i 



n— >oo 
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for all y, then dimsf >2 — a. 

Proof: See appendix A 

Note that part (a) of the theorem above is the generahzation of the statement that 
functions are locally Lipschitz (hence their graphs have dimension 1) to the case 

when the function has Holder type upper bound (hence their dimension is greater 

than one). 

Here the function is required to have the same critical order throughout the inter- 
val. We can weaken this condition slightly. Since we are dealing with box dimension 
which is finitely stable ||^, we can allow finite number of points having different crit- 
ical order so that we can divide the set in finite parts having same critical order in 
each part. 

We can also prove a partial converse of the above theorem. 

Theorem 6 Let f : [0,1] M be a continuous function. 

a) Suppose 

\fix)-fiy)\<c\x-yr, 

where c>0, 0<a<l and \x — y\ < S for some S > 0. Then 

y d\f{x) - f{y)) 

lim — = tor o < a 

^^y [d{x - y)]'^ 

for all y G (0, 1) 

b) Suppose that for each y G (0, 1) and for each 6 > there exists x such that 
\x — y\ < 6 and 

|/(x)-/(i/)|>cr, 
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where c > 0, 6 < 6q for some 5o > and < a < 1. Then there exists a sequence 
Xn y as n oo such that 



lim — -— -; = ±oo tor q> a 



for all y. 



Proof 

a) See appendix B 

b) Proof follows by the method used in the previous section to show that the 
fractional derivative of order greater than 2 — a of the Weierstrass function does not 
exist. 

These two theorems give an equivalence between Holder exponent and the critical 
order of fractional differentiability. 

.1 Appendix A: Proof of theorem 1 

(a) Without loss of generality assume = and /(O) = 0. We consider the case of 
q < a. 

As < g < 1 and /(O) = we can write 

c/-« d'if{x) 



fix) — 
^ ' dx-1 dxi 

jdy (0.17) 



T{q)Jo (x - y) i 
Now 

I difjy) I 
I dyi 

T{q)Jo ix^-i+^ 
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As, by hypothesis, for q < a, 



we have, for any e > 0, a 5 > such that \^--Pf^\ < e for all x < 5 



1/(^)1 < 



dxi 

e /"^ dy 



T{q)Jo 
e 



T{q+1 
As a result we have 

\f{x)\ < K\x\'' for x<S. 

Now this argument can be extended for general y simply by considering x — y instead 
of X and f{x) — f{y) instead of f{x). So finally we get for g < a 

\f{x)-f{y)\ < K\x-y\'^ ioi \x - y\ < 6, (0.18) 

for all y G (0, 1). Hence we have P 

dimB/(x) <2 — a. 

b)Now we consider the case q > a. If we have 

lim„^^^ = oo (0.19) 

3=71 ->0 dXn 

then for given Mi > and 5 > we can find positive integer such that \xn\ < S 
and d'^ f{xn)/dx'^ > Mi for all n > N. Therefore by eqn( p.l7| ) 

^^"^"^ - r{q)Jo (x„-i/)-«+i 
Ml , 
T{q + lf- 

68 



If we choose 5 = xn then we can say that there exists x < 5 such that 

f{x) > h6'^ (0.20) 

If we have 

Xn^O dXn 

then for given M2 > we can find positive integer N such that d'^f{xn)/dx'^ < —M2 
for all > A^. Therefore 

/(.„) < -"^r- iv 



T{q) Jo {xn - 
-M2 



^9 



r(g + i) - 

Again if we write 6 = x^, there exists x < 6 such that 

fix) < -M" (0.21) 
Therefore by ( p.20| ) and ( |0.21| )there exists x < 6 such that, for q > a, 

\f{x)\ > KS'^. 

Again for any y G (0, 1) there exists x such that for q > a and |a; — ?/| < S 

\f{x)-f{y)\ > k6^. 

Hence we have 

dimsfix) > 2 — a. 
69 



.2 Appendix B: Proof of the theorem 2a 



Assume that there exists a sequence x„ — ?/ as n — > cx) such that 
hm — -— --; = ±00 tor q < a , 



for some y. Then by arguments following the equation ( p. 19 ) it is a contradiction. 
Therefore 

hm — -— = const or tor g < a 

"^^y [d{x — y)\i 

Now if 

lim — 7^7-^ TT^-^ = const for q < a, 

"^^y [d[x — |/)J5 



then we can write 

dV{^) - f{y)) 



K + r]{x,y) 



[d{x - 

where K = const and ri{x, y) sufficiently fast as x y. Now taking e-derivative 
of both sides, for sufficiently small e we get 

d'-^'ifix) - fjy)) _ Kjx - y)~^ d^r]ix, y) 

[d{x - y)Y+^ r(l - e) + [d{x - y)Y °' ^ + ' " 

As X — s> y the right hand side of the above equation goes to infinity (term involving 
T] doesn't matter since 77 goes to sufficiently fast) which again is a contradiction. 
Hence the proof. 
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Chapter 5 



Pointwise Behavior of Irregular 
Functions 

5.1 Motivation 

The relevance of irregular functions in physical processes was discussed in Chapter ^ 
There are also cases such as intermittent bursts in velocity field of a turbulent fiuid, 
in chaotic signals etc. These and earlier examples demonstrate that the regularity of 
signals in physical systems may vary from point to point, and call for an elaborate 
study of pointwise regularity of irregular signals. Mult if r act al functions and measures 
are used to model these signals. 

There are several approaches to study pointwise behavior of functions. Recently 
wavelet transforms [|I|, |^ have been used for this purpose and have met with some 
success. In the following we demonstrate that LFD, as defined in chapter is a 
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tool that can be used to characterize irregular functions and point out its certain 
advantages over wavelet transforms. 



Recall that Holder exponent a{y) of a function f at y was defined in |2.4.1| as the 



largest exponent such that there exists a polynomial Pn{x) of order n that satisfies 

\f{x)-P^{x-y)\ = 0{\x-yn, (5.1) 

for X in the neighborhood of y. Various authors ^] have used this definition. It 
is clear from theorem |^ that LFDs provide a way to calculate Holder exponents and 
dimensions. It may be noted that since there is a clear change in behavior when order 
q of the derivative crosses critical order of the function it should be easy to determine 
the Holder exponent numerically. Earlier methods using autocorrelations for fractal 
signals ^ involve an additional step of estimating the autocorrelation itself. 

We begin our discussion on use of LFD to study pointwise behavior by consider- 
ing isolated and masked singularities. Then we present a short introduction on the 
multifractal functions and study their local properties using LFD. 



5.2 Isolated singularities and masked singularities 

Let us start 0, ^ with the case of isolated singularities. We choose a simple example 
f{x) = ax", 0<a;<l, x > 0, a constant. Critical order at x = gives the order 
of singularity at that point. 

For the cases where two or more singularities are superposed on each other at the 
same point, we can use LFD to detect a higher order singularity which may be masked 
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by one with lower order. As demonstrated below, we can estimate and subtract the 
contribution due to lower order singularity from the function and find out the critical 
order of the remaining function, in the following way. We first find the critical order 
of the given function at the point of the singularity. This itself is the order of the 
lower order singularity. Then we determine the value of the LFD at the critical order. 
Consider, for example, the function 

/(x) = az" + < a < p < 1, x > 0. (5.2) 

LFD of this function at x = of the order a is iD"/(0) = ar{a + 1). Using this 
estimate of stronger (lower order) singularity we now write 

G{x;a) = f{x)-f{0)-—^x" 

i (a + 1) 



which for the function / in eqn (|5.2| ) is 

d^Gjx'a) _ bTi/3 + 1) 
dxi T{P-q + l) 

Therefore the critical order of the function G, at x = 0, is /3. Notice that the estima- 
tion of the weaker (higher order) singularity was possible in the above calculation just 
because the LFD gave the coefficient of x'^/T{a + 1). This suggests that using LFD, 
one should be able to extract singularity spectrum masked by that of strong singu- 
larities. Hence one can gain more insight into the processes giving rise to irregular 
behavior. 

Comparison of two methods of studying pointwise behavior of functions, one us- 
ing wavelets and other using LFD, shows that characterization of Holder classes of 
functions using LFD is direct and involves fewer assumptions. 
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It has been shown [Q that by using wavelet transforms one can detect singular- 
ities masked by an nth order polynomial, by choosing analyzing wavelet with 
its first n moments vanishing. If one has to extend the wavelet method for 
the unmasking of weaker singularities from lower order ones, one would then 
require analyzing wavelets with fractional moments vanishing. If in addition 
there exists a polynomial masking both the singularities one will require the 
above condition along with the vanishing of first n moments. 

Using wavelets to characterize Holder class of functions with oscillating singu- 
larity, e.g., /(x) = x^sin{l/x^), x>0, 0<a<l and /5 > 0, will require two 
exponents |]16| . On the other hand LFD of the function, through the theorems ^ 



and directly gives its Holder exponent. 

• The class of functions to be analyzed using wavelet transforms is in general 
restricted . These restrictions essentially arise from the asymptotic properties 
of the wavelets used. On the other hand, with truly local nature of LFD one 
doesn't have to bother about behavior of functions outside our range of interest. 

Now we consider the functions where we have singularities of various orders at 
each and every point. 



5.3 Introduction to multifractal functions 



Extensive literature exists on the study of multifractal measures and their applica- 
tions IP, ^, |10|, |r^, [T^ . Their importance lies in the fact that such measures are natural 
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measures which are used in the analyses of many physical phenomena |T^, 14 1. How- 
ever there are cases in which the object one wants to understand is a function (e.g. 
a fractal or multifractal signal) rather than a set or a measure. For instance one 
would like to characterize the velocity field of fully developed turbulence. Moreover, 
to every measure we can assign a function (which is just its distribution function) but 
not vice versa. 



We saw in section [4.1|that the Weierstrass function is a fractal function, i.e., it has 



the same Holder exponent at every point. However, there exist more general functions 
which do not have same Holder exponent everywhere. These are the abovementioned 
multifractal functions. Such functions can be used to model various intermittent phe- 
nomena - spatial, temporal as well as spatio-temporal - arising in physical systems. 



Such an intermittent behavior is clearly seen in Fig. 5.1 which is a graph of a mul 



tifractal function. Recall that in section ^4.4|, while dealing with fractal functions, an 



equivalence between the Holder exponent and the critical order has been established. 
The same equivalence also holds for multifractal functions where we have different 
Holder exponents at different points. Selfsimilar multifractal functions have been 
constructed by Jaffard 0. Let us discuss one specific example of such a function. 
This function is a solution F[x) of the following functional equation 

d 

F{x) = J2 >^iFiS^\x)) + g{x), X real. (5.4) 

i=l 

where S'j's are the affine transformations of the kind Si{x) = fiiX + bi (with < 1 
and 6j's real), Aj's are some real numbers, and g is any sufficiently smooth function 
(it is assumed that g and its derivatives have fast decay) . For the sake of illustration 
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we choose the case d = 2, with /ii = /i2 = 1/3, bi = 0, 62 = 2/3, Ai = 3 A2 = 3 ^ 
(0 < a < /? < 1) and 

g{x) = sm(27rx) if x G [0, 1] 
= otherwise. 



The graph of F{x) with different values of parameters than above is plotted in Fig. 
Such functions are studied in detail in using wavelet transforms. There it was 
shown that the above functional equation ( |5.4| ), alongwith the parameters we have 
chosen above, has a unique solution F. Furthermore at every point F either has 
Holder exponents ranging from a to (3 or is as smooth a.s g. A sequence of points 
S'jj(O), Si^Sij^{0), ■ ■ ■ , Si^- ■ ■5'ji(0), ■ ■ ■, where takes values 1 or 2, tends to 
a point in [0, 1] (in fact to a point of a triadic Cantor set) and for the values of yUjS 
chosen above this correspondence between sequences and limits is one to one. The 
solution of the above functional equation is given by P 

00 2 

^(^) = E E K■■■>^^Ms^:■■■s^,\x)). (5.5) 

n=0 ii,---,in=l 

Note that with the above choice of parameters the inner sum in ( p. 51 ) reduces to a 
single term. Jaffard has shown that the local Holder exponent at y E [0, 1] is 

%)=liminfi^4^^M^l— (5.6) 

where {ii{y), ■ ■ -yiniy)} is a sequence of integers depending on a point y (i.e. Si-^^y) ■ 
■ ■ Si^(y) ^ y). It is clear that hmin = « and hmax = P- The function F at the points 
of a triadic cantor set has h{x) G [a,/?] and at other points it is smooth (as smooth 



as g). Benzi et. al. |15] have constructed multifractal functions which are random in 



nature unhke the above nonrandom functions. For still another approach see p5 



Several well-known pathological functions have been reanalyzed in |T8|, 0, ^ 
and found to have multifractal nature. Here we consider one classic example of 
multifractal function. 

oo 



1 



R{x) = ^ sinfvm^a; 



(5.7) 



This function was proposed by Riemann. It turns out that the regularity of this 



function varies strongly from point to point. Hardy and Littlewood ||2T[ proved that 
R{x) is not differentiable at Xq if Xq is irrational or if xq can not be written in the 
form 2p + l/2g+ 1 {p,q & N). In fact they showed that the Holder exponent at these 
points is less than or equal to 3/4. Gerver proved the differentiability of R{x) at 
points which can be written as 2p+ l/2g+ 1 {p,q G A^). At these points the exponent 



is 3/2. This function has also been studied in p3|, Jaffard IjT^l has recently shown 
that the dimension spectrum of R{x) is 



d{a) 



4a -2 if a G [|, f] 







if a 



— oo otherwise, 

where d{a) is the Hausdorff dimension of the set where the Holder exponent is a. 



5.4 Treatment of multifractal function 

We now proceed with the analysis of multifractal functions using LFD 0. Since LED 
gives the local behavior of the function, theorem ^ can also be applied to multifractal 
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functions. Jaffard's construction of a selfsimilar multifractal function has been studied 
in the last section. We consider the example ( |5.5| ) and study its LFD: 
d^jFjx) - F{y)) _ - 2 

d^igjS-J ■ ■ ■ Sr\x)) - gjSr' ■ ■ ■ S-\y))] 
[d{x - y)]" 

oo 2 
n=0 ii,---,j„=l 

dn9{S-:---S-\x))-giS-'---S-\y))] 

[d{S-: . . . S-\x - y))]'> ' ^'-'^ 

provided the RHS is uniformly bounded. Following the procedure described in sec- 
tion ^.2| the fractional derivative in the RHS can easily be seen to be uniformly 
bounded and series is convergent if g < mm{h{x),h{y)}. Further it vanishes in the 
limit X — *■ Therefore if g < h{y), lD'^F{y) = 0, as in the case of Weierstrass 
function, showing that h{y) is lower bound on critical order. 

The procedure of finding an upper bound is lengthy. It is relegated to the Ap- 
pendix below where it is shown that upper bound is also h{y). 

In this way an intricate analysis of finding out the lower bound on Holder exponent 
has been replaced by a calculation involving few steps. This calculation can easily be 
extended to cases incorporating more general functions g{x). 

Summarizing, the LFD enables one to calculate the local Holder exponent even 
for the case of multifractal functions. This fact, proved in theorems ^ and|^, has been 
demonstrated with concrete illustration. 
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.1 Appendix 

Our aim in this is appendix is to find an upper bound on the critical order and hence 
on the local Holder exponent for the function of equation ( |5.5| ). Our procedure will be 
similar to that of Jaffard the only difference being that we take fractional derivative 
instead of wavelet transform. We proceed as follows. The defining equation for F[x) 
when iterated times gives 

Af-l 2 

F{^) = E E K ■ ■ ■ KAsrJ ■ ■ ■ sr,\x)) 

+ j: A,---A,,F(5-;---5ri(a;)) a;G[0,l]. (0.10) 



ii .■••,«iv=l 



We now consider 



d^Fjx) - F{y)) _ f , , . 



[d{SrK . . Sr\x - y))] 



1 



2 

-Q 



n,-"i*jv=l 

d^[FiSr^' . . . Sr\x)) - ■ • ■ Sr\y))] 



. (0.11) 



[diS~^' ■ ■ ■ S-\x - y))]^ 

Let us denote the first term on RHS by A and the second term by 5. In the 
following y E (0, 1). Choose such that 3^(^+1) < \x — y\ < 3^^. Denote 
An{j/) = Ki{y) ■ ■ ■ ^in{y)- For the values of /ijS we have chosen /ij^ ■ ■ -/ij^ = 3~"'. Now 
since g is smooth \g{x) — g{y)\ < C\x — y\. 

d''[g{Sr:---Sr\x))~giSrJ...Sr\y))] ^ x - y ^,„^ ^^^^^^ 

[d{Sr„^ ■ ■ ■ S^^\x - y))]i - 



Prom the way we have chosen /ijS and |x — ?/| this term is bounded by C3"^^ ''^3 ^'•^ 
Therefore the first term in equation ( |0.11| ) above is bounded by 



N-l 

= C3-^(i-'?) 5] A„(,)3" 

n=0 

< C3-^(i-'^)A(^_i)(,)3^-i(l + 



A 



A 



3A(jv_i)(2^) 32A(7v_i)(y) 



+ ■ ■ 



(0.13) 



The quantity in the bracket is bounded. Therefore 



N{y) 



(0.14) 



Now we consider the second term in the RHS of equation ( |0.11| ) and find a lower 
bound on that term. We assume that F{x) is not Lipschitz at y (for otherwise LHS 
of equation (|0.11| ) will be zero in the limit x ^ y and the case is uninteresting). 
Therefore there exists a sequence of points x„ approaching y such that 



\F{xn) - F{y)\ > c\xn - y\ 



(0.15) 



The function is not Lipschitz at y implies that it not Lipschitz at S^^^ ■ ■ ■ Si^{y). 
Therefore we get (one may recall that only one of the 2" terms in the summation 
SL-,iiv=i is nonzero) 



B > cA^(,)3^^3^(i-'^)3-^(i-'^) 



(0.16) 



Therefore there exists a sequence x„ such that 
^df^iF^x^) - F{y)) 



[d{xn - y)Y 



— cAAr(y)3^''| < C\N{y)'i^'^ 



(0.17) 
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Since (p.l5|) is valid for every c for large enough 



n 



d%F{xn) - F{y)) 



where C = c — C Now from the expression of h{y) it is clear that when q > h{y) 
LFD does not exist and the critical order is bounded from above also by h{y). 
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Chapter 6 



Local Fractional Differential 



Equations 



6.1 Motivation and Introduction 

The most important equations of physical sciences are differential or integral equa- 
tions, two themes emerging from the seventeenth-century calculus. The derivatives of 
integer order are coefficients of integer powers in the Taylor expansion. In particular 
the first order derivative models a local linear behavior and the existence of the first 
order derivative of the function implies Lipschitz continuity (i.e., the Holder exponent 
> 1 as defined in |2.4.1j ) of the function. Therefore it is not possible to have Holder 



continuous functions as solutions to ordinary differential equations. Such functions 
arise while dealing with phenomena involving fractals. Hence it becomes necessary 
to develop calculus which would go beyond ordinary calculus and incorporate Holder 
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continues functions too. 

The local fractional derivatives that we have studied in previous chapters preserve 
the local nature of the derivative operator. Moreover they characterize the Holder 
exponent of an irregular function. This suggests these LFDs may provide a much 
needed tool for doing calculus for fractal space-time. Hence it becomes pertinent to 
ask the questions such as: what is the inverse of the LFD operator (if it exists), or can 
one write and solve equations involving LFDs? Also, what meaning and applications 
these local fractional differential equations will potentially have? These questions 
open a big area. We'll try to answer some of these questions in the following sections. 

6.2 Local fractional differential equation 

In this section we introduce a new type of equation which involve LFDs. We call them 
local fractional differential equations (LFDE). These are equations involving LFD of 
a dependent variable with respect to independent variables, as in the case of ordinary 
differential equations. In order to understand the meaning of such equations let us 
consider a simple example of the LFDE: 



In the above equation g[x) is a known function and we have to find f{x) such that 
this equation is satisfied. The issue of the conditions guaranteeing solutions of such 
an equation is too involved to be addressed here in its full generality. Therefore we 
restrict ourselves to specific cases. 




(6.1) 
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We begin by noting that the equation 

E>lf{x) = Const, (6.2) 



does not have a finite solution when < g < 1. Interestingly, the solutions to ( |6.1| ) 
can exist, when g{x) has a fractal support. For instance, when g{x) = xc{x), the 
membership function of a cantor set C (i.e. g{x) = 1 if x is in C and g{x) = 
otherwise), the solution with initial condition /(O) =0 exists if g = a{= dim^/C). 
Explicitly, generalizing the Riemann integration procedure. 



i=0 



r(a + l) ^ r(a + l)^ 



where subdivision points of the interval [xq = 0, X]\[ — X ] and is a fiag 

function which takes value 1 if the interval [xj, Xj+i] contains a point of the set C and 
otherwise. Now if we devide the interval [0, x] into equal subintervals and denote 
A = Aj = Xj+i — Xi then we have Pc{x) = A^X^-^c- But J2 order of 

iV~°. Therefore Pc{x) satisfies the bounds ax" < Pc{x) < where a and 5 are 
suitable positive constants. Note that Pc{x) is a Lebesgue-Cantor (staircase) function 
as shown in the Fig. |6.1[ Since the function Xc'{x) is zero almost everywhere the 
function Pc{x) is constant almost everywhere. As is clear from the equation ( |6.3| ), it 
rises only at points where Xc{x) is non-zero. In general, the algorithm of the equation 
( |6.3D is expected to work only for the sets C for which dim^C = dimj^C (in fact in 
this case only N°' terms in the summation are nonzero). It can now be seen that if 
g{x) = Const then the sum in equation ( |6.3| ) does not exist. Therefore the equation 
ID%f{x) = Const does not have a solution as stated earlier. Though this algorithm 



0.2 0.4 0.6 0.8 1 

Figure 6.1: A Lebesgue-Cantor staircase function 

possibly can be extended for more general sets we choose to restrict ourselves to the 
sets satisfying the above condition. 

It may be noted that the above LFDE is different from usual fractional differ- 
ential equations. Conventional fractional differential equations, i.e., equations in- 
volving Riemann-Liouville fractional derivatives, have been studied in the literature. 
These equations have found applications at various places, e.g., solution of Bessel's 



equation, diffusion on curved surface, etc [O. In fact the equations appearing in 



0, ^ ^ 1^, 1^, ^, |TT| are of this type. Since the Riemann-Liouville fractional deriva- 
tive can be written as a integral of some function these equations are essentially 
integro-differential equations. On the other hand the equation we have proposed is 
intrinsically local. As mentioned in chapter 0, LED E)'^f{y) gives the coefficient of 
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the power {x — yY in the fractional Taylor expansion and hence local fractional dif- 
ferential equations are equations satisfied by these coefficients. We feel that this is a 
more natural generalization of ordinary differential equations. This is not only due to 
the fact that the local nature is preserved but also because the ordinary differential 
equations can be interpreted in similar manner, i.e., they are equations among the 
coefficients of integer powers (x — y)" in the expansion of a function f{x) at y. 

The equation ( |6.1| ), along with its solution (|6.3|) for the special choice of g{x), has 
an interesting consequence. When q = a = 1, i.e. the set C has dimension 1, we get a 
equation /'(x) = Xc{x). The solution f{x) of this equation gives the length of the set 
[0, x] n C. In the similar manner, the solution of the equation ( |6.1| ) with g{x) = Xc{x) 
would yield the fractal measure of the set [0, x] f] C. This indicates that the operation 
of the inversion of the LFD, a 'local fractional integral' (not to be confused with 
fractional integral), could yield a definition of the fractal measure. It is interesting 
to note that this definition, unlike other definitions of fractal measures, would have 
its origin in calculus. This definition is equivalent, atleast in the case of simple sets 
considered here, to the Hausdorff measure except for the normalization constant. The 
normalization constant l/r(a + 1) in equation ( |6.3D comes from calculus whereas the 
one in Hausdorff measure (equation (|2.34| )) comes from the geometry of the covering 
elements. Of course, all these observations regarding the equivalence between two 
approaches, other than the normalization constant, remains to be rigorously proved 
for general sets. But this interesting observation makes it plausible to answer the 
questions raised in the chapter ^ about the relation between 'a-fold' integral and the 
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a-dimensional measure. 



6.3 Derivation of local fractional Fokker-Planck equa- 
tion 

The basic idea of this section is to utihze the fractional Taylor expansions, developed 
in chapter ^ in the Chapman-Kolmogorov condition and obtain analogs of Fokker- 
Planck-Kolmogorov (FPK) equation in the LFD calculus formalism. We briefly re- 
call the usual procedure of deriving the FPK equation Let W{x,t) denote the 



probability density for a random variable X taking value x at time t. Then the 
Chapman-Kolmogorov equation is 

W{x, t + T)= [ P{x, t + r|x', t)W{x', t)dx\ (6.4) 



where P(a;i, ti|a;2, ^2) denotes the transition probability from X2 at time ^2 to Xi at 
time ti and r > 0. In order to obtain the usual FPK equation we first expand 
the integrand in the equation ( |6.4D in a Taylor series with respect to x — x' . Further 
we assume that the moments, 

Mn{x, t,T) = I dy{y - xyP{y, t + t\x, t), (6.5) 

can also be expanded into Taylor series with respect to r as 

^"^""'/'^^ = yl"(a;, t)r + 0{t^). (6.6) 
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Then various powers of r are equated by expanding LHS with respect to r. As a 
result one arrives at the FPK equation 

^^^ = C{x,t)W{x,t), (6.7) 

where the Fokker-Planck operator C{x^t) is 

C{x,t) = - A\x,t) + —A\x,t). (6.8) 

ox OX"^ 

There are a number of hmitations of this approach because of various assumptions 
going into its derivation. For instance, as noted in |[T^, probabihty distributions whose 



second moment does not exist are not described by FPK equation even though such 
distributions may satisfy original Chapman-Kolmogorov equation. Also, as empha- 



sized in 1 18], the differentiability assumption going into the Taylor expansion of the 



integrand of equation (|6.4| ) may also break down in various situations. For instance, 
the transitional probability density may not be differentiable at x = x'; in which case 
the derivation of FPK equation itself will break down. Another situation arises when 
we have a fractal (nondifferentiable) function as the initial probability density. In 
such a case even the usual Fokker-Planck operator can not be operated on the initial 
density. 

It is thus of interest to broaden the class of differential equations one can derive 
starting from the Chapman-Kolmogorov equation, and study various processes de- 
scribed by them. In this section we investigate the possibility of dealing with those 
probability densities which need not satisfy usual differentiability criteria. We follow 



the standard method to derive the Fokker-Planck equation from equation (|6.4| ) out- 
lined above with the modification that we now expand the integrand using fractional 
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Taylor expansion ( p.l4| ) instead of ordinary Taylor expansion. Thus, if A = x — x', 
we have to consider two cases A > and A < separately. When A > 

P{x,t + T\x',t)Wix',t) = P{x - A + A,t + T\x - A,t)W{x- A,t) (6.9) 

^-—A^E>^[P{x + A,t + T\x,t)W{x,t)], (6.10) 



r(/9 + 1) 



and when A < 

P{x,t + T\x',t)W{x',t) = P{x - A + A,t + T\x - A,t)W{x - A,t) (6.11) 

^ 1 d 
^0^(^ + 1) d{-x) 

+^;^^^^(-A)''iD^[P(x + A, t + r|x, t)W^(x, t)],(6.12) 

where N < (3 < N + 1. This is the only place where our derivation differs from that of 
conventional Fokker-Plank equation. (In the conventional derivation one expands the 
integrand in equation ( |6.4| ) in Taylor series assuming implicitly that the probability 
densities involved are analytic.) Substituting this into equation ( |6.4| ) we get 

^1 Br 

W{x,t + T) = W{x,t) + Y,——-(-^r dx'A''P{x + A,t + T\x,t)W{x,t) 

^i^{n + l) d{-x) J 

^ r(/+l) ^^-^/.°° - x)''P{y, t + r\x, t)W{x, t)] 

+ — — - — — -iDf+[ / dy{x — y)^P{y, t + t\x, t)W{x, t)] +Remainder 
i yfj ~l~ Ij J —oo 

where E)^ is a partial LFD w.r.t x. Now if < a < 1 

W(x, t + r) - W(x, t) = + Remainder. 

T(a + 1 
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where iDj is partial LFD w.r.t. t. In general a and /? may depend on x and t. But 
we assume that a and /3 are constants. Therefore we get 

r(« + i) " tid{-x)^ ^ r(n + i) 

Mt(x,t,T) , 

+iD.M ^g";';^V (x,t)i (6.13) 

where 

/•OO 

M+{x,t,r)^ dy{y-xyP{y,t + r\x,t) a > 0, (6.14) 

Jx 

M-{x,t,T)=r dy{x -yyP{y,t + T\x,t) a>0 (6.15) 

J—oo 



and 



M„(x, t, r) = M+(a:, t, r) + M-(x, t, r) (6.16) 

are transitional moments. The limit r — > gives us the local fractional Fokker- 
Planck-Kolmogorov (FFPK) equation 

E)fW(x, t) = jC^(x, t)W(x, t) (6.17) 

where the operator is given by 

AT o 

Ci{x,t) = ^^-^)-Al{x,t) + ]DtAL{x,t) + ]Dl^Ai^{x,t) (6.18) 

n=l 

where 

^2,(..0.1im^fc^n^^ (6.19) 
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and 

A^^{x,t) = Al+{x,t)+Ai_{x,t). (6.20) 

Here corresponding A^s are assumed to exist. We would like to point out that the 
equation (|6.17|) is analogous to truncated Kramers-Moyal expansion. Two rather 



important special cases are, < /? < 1 and 1 < (3 < 2. In the former case we get the 
operator 

C^,{x,t) = ]DtAi.{x,t)+lD^,+Ai+{x,t), 0</3<l, 

and in the latter case we get 

Ci{x,t) = -^Al{x,t)+lDtAL{x,t)+lDl^Ai^{x,t), l<(3<2. 

This operator can be identified as generalizations of the Fokker-Planck operator (in- 
volving one space variable). It is clear that when a = 1 and /5 = 2 we get back the 
usual Fokker-Planck operator in ( |6.8|) . 

Returning back to expression (|6.19|) it is clear that the small time behavior of 



different transitional moments decide the order of the derivative with respect to time 
(in order to demonstrate this point we consider the example of a Levy process below). 
On the other hand, small distance behavior of transitional probability or the differ- 
entiability property of the initial probability density would dictate the order of space 
derivative. Depending on the actual values of a and (3 as well as their interrelation 
the above local FFPK equation will describe different processes. 

Equations which give rise to evolution-semigroup are of interest in physics. The 
equation (|6.171 ) corresponds to a semigroup if a = 1. One can then write down a 
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formal solution of the above equation in this follows. In the time independent 

case we have 

iy(a;,t) = e^i(^)*iy(x,0) (6.21) 
and when the operator depends on time we have 

W{x, t) =T e/o 0) (6.22) 

where £f is an operator in equation (|6.17D and T is the time ordering operator. 



6.4 Examples 

Levy processes 

The symmetric stable Levy process of index /x, < < 2, is defined by 

/oo 
dke-^\^\"'e'^\ (6.23) 
-oo 

For this process all the moments, 

/•oo /'OO 

M+(t) = / dxx^ / dA;e-''l^l"*e^'=^ (6.24) 
of order 7 > /i diverge (see [0). We have 

/•oo roo 

M+{\t) = / dxx'' / dke-^\^\'^^'e'^'' (6.25) 

' JO J-oo 

= A'^/"M+(t). (6.26) 

A similar relation holds for M~ . As a result the moments scale as M^{Xt) = X^^^M^it) 
and the assumption in equation ( |6.6| ) is clearly not valid. If we use our local FFPK 
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then we get 

lD]^^W{x,t) = lD24A^^^^_{x,t)W{x,t)]+lD24A]^^4x,t)W{x,t)]. (6.27) 

Since the process is symmetric the first derivative does not appear. The order of the 
time derivative depends on that of space derivative but it is always less than one. 
There is only one free parameter 7 which is restricted to be in the range < 7 < 
fi. In this case the value of 7 will be decided by the differentiability class of the 
initial distribution function. When /i = 2 we get a Gaussian process whose second 
moment is finite; hence we can have 7 = 2 in the above equation and get back 
the usual Fokker-Planck equation. Equation ( |6.27| ) forms one example where we get 



nontrivial (fractional) values for the orders of the derivatives and usual derivation of 
FPK equation breaks down. 

A fractal time Gaussian process 

As our next example we consider the transition probability 



P{x,t + T\x',t) = , - e^Pc(t.r) (6.28) 

/7rAPc(t,r) 



= 6{x-x') if APc{t,T) = 0, (6.29) 

where APc{t, r) = Pc{t + T) —Pc{t)). This transition probability describes a nonsta- 
tionary process which corresponds to transitions occurring only at times which lie on 
a fractal set. Such a transition probability can be used to model phenomena such as 
diffusion in the presence of traps. The second moment is given, from equation (|6.16| ), 
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by 

^ ' ^ 2 2r(a + l) 

= yXcW (6.30) 

This gives us the following local fractional Fokker-Planck equation (in this case an 
analog of a diffusion equation). 

E>tWix,t) = ^^i^plxcit)-^Wix,t) (6.31) 

We note that even though the variable t is taking all real positive values the ac- 
tual evolution takes place only for values of t in the fractal set C. The solution of 
equation (|6.31| ) can easily be obtained as 



W{x,t) = Pt.toW{x,to) (6.32) 

where 

iV-l 1 Q2 

The above product converges because except for terms for which = 1 (which are 
of order A^") all others take value 1. It is clear that for to < t' < t 

Wix,t) = Pt^t,Pt,_t,W{xM) (6.34) 

and Pt gives rise to a semigroup evolution. Using equation ( |6.3| ) it can be easily seen 
that 

W{x, t) = e^^W{x, to = 0). (6.35) 
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Now choosing the initial distribution W{x, 0) = 6{x) and using the Fourier represen- 
tation of delta function, we get the solution 



W{x,t) = (6.36) 



Consistency of the equation ( |6.3(: ) can easily be checked by directly substituting this 



in Chapman- Kolmogorov equation (|6.4| ). We note that this solution satisfies the 
bounds 

^ -.e^ <W{x,t) < ^^=e^ (6.37) 



for some < a < 6 < oo. This is a model solution of a subdiffusive behavior. It 
is clear that when a = 1 we get back the typical solution of the ordinary diffusion 
equation, which is (vrt)"-*^/^ exp(— x^/t). 

Anomalous diffusive behavior arises at numerous places in physics and therefore 
obtaining such solutions from a kinetic equation is very important. Specifically such 
transition probabilities (equation (|6.28|)) may arise in disordered systems when there is 



a diffusion in presence traps, diffusion across cantori in phase space of the Hamiltonian 
systems or in phase space of weak chaotic Hamiltonian systems. 

This work may have interesting consequences in other fields also. The LFDEs 
that we have considered above suggest a generalization of Dynamical systems to 
orders less than one. The solution that we have found is also a solution of the 
Chapman- Kolmogorov equation defining Markov process. Therefore our results may 
have relevance in probability theory. We would like to suggest that we have been able 
to incorporate the phenomenon taking place in fractal time into an equation and find 
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its exact solution. This marks a step towards development of calculus for phenomena 
occurring in fractal space-time. 
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Chapter 7 



Concluding Remarks 



• The thesis introduced and discussed a new notion called local fractional deriva- 
tive of a function of one variable. This is a modification of conventional frac- 
tional derivatives. A quantity, critical order of a function at a given point, has 
been defined as the maximum order such that all the local fractional derivatives 
of orders below it exist. 

• A fractional Taylor expansion has been constructed which contains these LFDs 
as coefficients of fractional power law. This gives a geometrical interpretation 
to LFD. 

• We have established that (for sufficiently large A ) the critical order of the 
Weierstrass' nowhere differentiable function is related to the box dimension of 
its graph. If the dimension of the graph of such a function is 1 + 7, the critical 
order is 1 — 7, i.e., the function is differentiable upto order 1 — 7. When 7 
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approaches unity the function becomes more and more irregular and local frac- 
tional differentiability is lost accordingly. Thus there is a direct quantitative 
connection between the dimension of the graph and the fractional differentiabil- 
ity property of the function. A consequence of our result is that a classification 
of continuous paths (e.g. fractional Brownian paths) or functions according to 
local fractional differentiability properties is also a classification according to 
dimensions (or Holder exponents). 

• Also the Levy index of a Levy flight on one dimensional lattice is identified 
as the critical order of the characteristic function of the Levy walk. More 
generally, the Levy index of the distribution is identified as the critical order of 
its characteristic function at the origin. 

• We have argued and demonstrated that LFDs are useful for studying isolated 
singularities and singularities masked by the stronger singularity (not just by 
regular behavior). We have further shown that the pointwise behavior of irreg- 
ular, fractal or multifractal, functions can be studied using LFDs. 

• We have also demonstrated that it is possible to carry out the same theme even 
in the multidimensional case. In particular, the Holder exponents in any direc- 
tion are related to the critical order of the corresponding directional-LFD. We 
feel that, whereas a one dimensional result is useful in studying irregular signals, 
the results here may have utility in image processing where one can characterize 
and classify singularities in the image data. We note that it is possible to write 
a multivariable fractional Taylor series which can be used for approximations 
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and modeling of multivariable multiscaling functions. Further these considera- 
tions suggest a way for formulating fractional differential geometry for fractal 
surfaces. 

• We have introduced new kind of differential equations called local fractional 
differential equations. They are equations in terms of LFDs. Our equations are 
fundamentally different from any of the equations proposed previously, since 
they involve LFDs. They are local and more natural generalization of ordinary 
differential equations. LFDEs deserve a thorough study in their own right. 

• We have derived a generalization of FPK equation which involves the local 
fractional derivatives. Since the present analog of FPK equation is derived 
from first principles, we feel that our equations will have general applicability 
in the field of physics. We expect them to be of value in the studies of anomalous 
diffusion, disordered phenomena, chaotic Hamiltonian systems, etc. 
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